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j_i ' Abstract. 

Qh' In the current paper, we further develop the model for the migration 

I of planets introduced in Del Popolo et al. (2001) and extended to time- 

dependent planetesimal accretion discs in Del Popolo & Ek§i (2002). We 
use a method developed by Stepinski & Valageas (1996, 1997), that is 
fSI ■ able to simultaneously follow the evolution of gas and solid particles for 

I up to lO^yr. The disc model is coupled to the migration model introduced 

in Del Popolo et al. (2001) in order to obtain the migration rate of the 
planet in the planetesimal disc. We find that in the case of discs having 

■ total mass of 10~^ — O.IMq, and 0.1 < a < 0.0001, planets can migrate 
\ inward a large distance while if M < 10~^Mq the planets remain almost 

■ in their initial position for 0.1 < a < 0.01 and only in the case a < 0.001 
\ the planets move to a minimum value of orbital radius of ~ 2AU. The 

■ model gives a good description of the observed distribution of planets in 
the period range 0-20 days. 
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' 1. Introduction 

H . 

. 5t , In order to study the formation of planetary systems it is necessary to study 

the global evolution of solid material which constitutes, together with gas, the 
protoplanetary discs. In the following, we introduce a time-dependent accretion 
disc model that shall be used in the next sections to study planets migration. 

It is usually assumed that the surface density in planetesimals is propor- 
tional to that of gas (Lin & Papaloizou (1980, page 47); Murray et al. (1998)) 
and that the spatial structure and the time evolution of the surface density, 
T,s{R,t), is such that 'Es{R,t) oc S(i?, t). Although this approximation has been 
and is widely used in literature, a better way to proceed is to directly calculate 
the evolution of solids in the disc. We know that the evolution of the surface den- 
sity of gas is described by a diffusive- type equation, while that of solid particles 
is an advection-diffusion equation: 
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dt r dr 
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where fg = gr, the Schmidt number, Sc, is given by: 



sc = (i + j^k4)yi + |F2 (2) 

where v is the relative velocity between a particle and the gas, Vt the turbulent 
velocity, is the Keplerian angular velocity, tg the so called stopping time. 

Stepinski &; Valageas (1996, 1997) developed a method that, using a series 
of simplifying assumptions, is able to simultaneously follow the evolution of gas 
and solid particles due to gas-solid-coupling, coagulation, sedimentation, and 
evaporation/condensation for up to lO^yr. 

The equation of gas evolution is solved by means of an implicit scheme and 
the evolution of gas is computed independently from the evolution of particles. 
At every time step the quantities needed for evaluating the change in the mass 
distribution of solids are calculated and the change itself is computed from Eq. 
(1) using the operator splitting method. In such a method the advective term 
in Eq. (1) is treated by the numerical method of characteristics, whereas an 
implicit scheme is applied to the diffusion term. The obtained distribution of 
solid material is then modifed because of the existence of the evaporation radius, 
and the mass distribu- tion of the vapor is calculated using the implicit scheme. 
Finally, the new particle size distribution is calculated be- fore proceeding to 
the next time step. 

We suppose that a single planet moves in a planetesimal disc under the 
influence of the gravitational force of the Sun. The equation of motion of the 
planet can be written as: 

= Fo + R (3) 

(Melita &: Woolfson 1996), where the term Fq represents the force per unit mass 
from the Sun, while R is the dissipative force (the dynamical friction term-see 
Melita & Woolfson 1996). 

In order to take into account dynamical friction, we need a suitable formula 
for a disc-like structure such as the protoplanetary disc. 

We assume that the matter-distribution is disc-shaped, then we have that: 

R = — A;||i;i||e|| — k±Vi±e± (4) 

where ey and ej^ are two versors parallel and perpendicular to the disc plane 
and fc|| and k± are given in Del Popolo & Ek§i (2002). 

2. Results 

Our model starts with a fully formed gaseous giant planet of IMj at 5.2 AU. 
According to several evidences showing that the disc lifetimes range from 10^ 
yr to 10^ yr (Strom ct al. 1993; Rudcn &: Pollack 1991), we assume that the 
disc has a nominal effective lifetime of 10® years (Zuckerman et al. 1995). We 
integrated the model introduced in the previous section for several values of disc 
masses: Md = 0.1, 0.01, 0.001, 0.0001 Mq, and several values of a. 
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The results of the disc model are plotted in Fig.l. Fig. la shows the 
evolution of the density of the gas component of the disc in the case M = O.IMq 
and a = 0.1 (see figure). As shown in Fig. lb, in agreement with Stepinski &: 
Valageas (1997) and Kornct et al. (2001), the most important result of the 
low-mass models calculation (discs with M < O.IMq, and angular momentum 
~ lO^^gcm^/s) is that such models lead to the survival of solid material, as 
can be seen from the evolution of the particle size, or by the emergence of the 
converged, nonvanishing surface density distribution of solids. 

In the successive figures (2a-2d), we plot the evolution of semi-major axis 
of the planet (see figure). Summarizing figures 2a-2d, according to the final 
distribution of planets distances, the present model predicts that, unless the disc 
mass is very small M ~ O.OOOIMq, planets tend to move close to the central 
star to distances of the order of 0.03AU. It is also evident that it is possible to 
find a planet at any distance from their locations of formation and very small 
distances from the parent star for peculiar values of the parameters a and M. 

In Fig. 3a, in order to show the predictions of the model for what concerns 
the distribution of planets in the inner part of the disc, we plot the fraction of 
planets in the orbital period range 0-20 days calculated using the model of this 
paper. Fig. 3b, represents the same distribution obtained with the data given 
in www.cxoplancts.org (sec also Kuchner & Lecar 2002). 

According to Kuchner &; Lecar 2002, disc temperature determines the or- 
bital radii of the innermost surviving planets, similarly to our model. As streseed 
in DPI, DP2, the model has not the drawback of Murray et al. 1998 model, 
namely that of requiring a too large disc mass for migration and at the same time 
has the advantage of Murray et al. 1998 model of having an intrinsic natural 
mechanism that provides halting of migration. 

References 

Del Popolo, A., Gambera, M., E. N. Ercan, 2001, MNRAS, 325, 1402 

Del Popolo, A., Ek§i, Y., E. 2002, MNRAS, 332, 485 

Kornet, K., Stepinski, T. F., Rzyczka, M, 2001, A&A 378, 180 

Kuchner, M.J., Lecar, M., 2002, astro-ph/0206232 

Lin D.N.C, Papaloizou J., 1980, MNRAS, 191 

Melita M.D., Woolfson M.M., 1996, MNRAS 280, 854 

Murray, N., Hansen, B., Holman, M., Tremaine, S., 1998, Science 279, 69 

Ruden, S.P, Pollack, J.B., 1991, ApJ 375, 740 

Stepinski, T.F., Valageas, P., 1996, A&A 309, 301 

Stepinski, T.F., Valageas, P., 1997, A&A 319, 1007 

Zuckerman, B., Forveille, & Kastner, J.H., 1995, Nature 373, 494 



4 



A. Del Popolo 




100 1000 10' 100 1000 

(2c) (2d) 




Figure 1. (a) Evolution of gas for a disc with Md = O.IMq and 
a = 0.1 at t = lO^yrs (dashed Hne), t = lO^yrs (dotted Hne) and 
t = lO^yrs (soUd hne). (b) Same as Fig. la but for sohds. 

Figure 2. (a) The evolution of the a(t) of a Jupiter-mass planet, 
M = IMj in a planetesimal disc for a = 0.1 and several values of 
Md, 0.1 (sohd line), 0.01 (dotted line), 0.001 M© (short-dashed line) 
and 0.0001 (long-dashed line). (b,c, d) Same as Fig. (6a) but with 
a = 0.01,0.001,0.0001 

Figure 3. (a) Fraction of planets having orbital periods in the range 



